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^ ■ Abstract 

W. 

We consider the transverse-momentum (qx) distribution of Drell- Yan lepton pairs 
produced in hadron collisions. At small values of qx, we resum the logarithmically- 
enhanced perturbative QCD contributions up to next-to-next-to-leading logarithmic 
accuracy. At intermediate and large values of g*r, we consistently combine resummation 
with the known next-to-leading order perturbative result. All perturbative terms up to 
order a\ are included in our computation which, after integration over g^, reproduces 
the known next-to-next-to-leading order result for the Drell- Yan total cross section. 
We show and discuss the reduction in the scale dependence of the results with respect 
to lower-order calculations, estimating the corresponding perturbative uncertainty. We 
present a preliminary comparison with Tevatron Run II data. 
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1 Introduction 

The hadropro duct ion of the vector bosons W and Z/j*, also known as the Drell-Yan (DY) 
process [I], is the process to which parton model ideas (previously developed for deep inelastic 
lepton-hadron scattering) were first applied in the context of hard-scattering processes in hadron- 
hadron collisions. 

At high-energy hadron colliders, such as the Tevatron and the LHC, vector bosons are produced 
with large rates and with relatively-simple experimental signatures. The vector boson production 
process is thus relevant for various reasons. It is important for detector calibration; it provides 
us with strong tests of perturbative QCD and, in particular, it gives stringent information on the 
parton densities of the colliding hadrons; it represents an important background for new-physics 
searches. Owing to these reasons, it is essential to have accurate theoretical predictions for vector 
boson production cross sections and related kinematical distributions. 

These predictions are based on perturbative QCD and are obtained as power series expansions 
in the strong coupling as- The total cross section [2] and the rapidity distribution of the vector 
boson [3] are known up to the next-to-next-to-leading order (NNLO) in QCD perturbation theory. 
Fully exclusive NNLO calculations, including the leptonic decay of the vector boson, are also 
available jU El E]. Electroweak corrections up to 0(a) have been computed both for W [7j and 
Z/j* production [8]. 

In this paper we consider the transverse-momentum (q?) spectrum of the vector boson. The qx 
spectra of the W and Z bosons are particularly important since the uncertainties in their shape 
directly affect the measurement of the W mass. In the large-g-r region (q? ~ my), where the 
transverse momentum is of the order of the vector boson mass my, perturbative QCD calculations 
based on the truncation of the perturbative series at a fixed order in a$ are theoretically justified. 
In this region, the QCD radiative corrections are known up to the next-to-leading order (NLO) 
[9j [TOj [TI]. Nonetheless the bulk of the vector boson events is produced in the small-g-r region 
(<7t <C m 'v)i where the convergence of the fixed-order expansion is spoiled by the presence of large 
logarithmic terms, ag \n m (my/q^). To obtain reliable predictions, these logarithmically-enhanced 
terms have to be systematically resummed to all perturbative orders p2]-[20]. The resummed 
and fixed-order calculations at small and large values of qx can then be consistently matched at 
intermediate values of qx, to obtain QCD predictions for the entire range of transverse momenta. 

We use the transverse-momentum resummation formalism proposed in Refs. [20| [2Tj . The 
formalism is valid for a generic process in which a high-mass system of non strongly-interacting 
particles is produced in hadron-hadron collisions. The method has so far been applied to the 
production of the Standard Model (SM) Higgs boson [221 1211 [23] , single vector bosons [21], WW 
[25] and ZZ [26J pairs, slepton pairs [22], and DY lepton pairs in polarized collisions [2H]. The 
study of Ref. [21] is mainly based on next-to-leading logarithmic (NLL) resummation at small qr 
and on the leading-order (LO) calculation at large q?. In this paper we extend the analysis and 
the results of Ref. [21], combining the most advanced perturbative information that is available 
at present: next-to-next-to-leading logarithmic (NNLL) resummation at small qx and the NLO 
calculation at large qr- Other phenomenological studies of the vector boson qr distribution, which 
combine resummed and fixed-order perturbative results at various levels of theoretical accuracy, 
can be found in Refs. 



The paper is organized as follows. In Sect. H]we briefly review the resummation formalism of 
Refs. [20j EI] and its application to vector boson production. In Sect. [3] we present numerical 
results for Z/j* production, and we comment on their comparison with the Tevatron Run II data 
[30l l3T] . We also study the scale dependence of our results to the purpose of estimating the 
corresponding perturbative uncertainty. In Sect. H]we summarize our results. 



2 Transverse- momentum resummation 



We briefly recall some of the main points of the transverse-momentum resummation formalism of 
Refs. [201 EI]- Here we consider the specific case of DY lepton pair production, i.e. the production 
of a vector boson V (V — W + , W~ , Z/j*) that subsequently decays in a lepton pair. 

The inclusive hard-scattering process is 

h 1 (p 1 )+h 2 (p 2 ) -»- V(M,q T )+X -> h + h + X, (1) 

where h% and h 2 are the colliding hadrons with momenta p\ and p 2 , V is the vector boson (which 
decays in the lepton pair Z l5 1 2 ) with invariant mass M and transverse momentum q?, and X is an 
arbitrary and undetected final state. 

According to the QCD factorization theorem the q F differential cross section dayfdq\ can be 
written as 

-tt(Qt,M,s) = Y] / dx 1 / dx 2 f a/hl (x 1 ,fxl)f b/h2 (x 2 ,fj 2 F )—^(q T ,MJ;a s (fxj i ),fx 2 R ,fx 2 F ) , 
aq T ab Jo Jo ac lT 

(2) 
where f a /h(%, Hf) ( a = Qi Qj d) are ^ e parton densities of the colliding hadron h at the factorization 
scale fj,p, da^ b /dq T are the perturbative QCD partonic cross sections, s (s — Xix 2 s) is the square 
of the hadronic (partonic) centre-of-mass energy, and ixr is the renormalization scale. 

In the region where q? ~ M (in practice, we always consider the case in which M is close to 
the mass my of the vector boson), the QCD perturbative series is controlled by a small expansion 
parameter, as(M), and fixed-order calculations are theoretically justified. In this region, the QCD 
radiative corrections are known up to next-to- leading order (NLO) [9]. 

In the small-gT region (q? <C M), the convergence of the fixed-order perturbative expansion is 
spoiled by the presence of powers of large logarithmic terms, a^ In 171 (M 2 /q F ). To obtain reliable 
predictions these terms have to be resummed to all orders. 



We perform the resummation at the level of the partonic cross section, which is decomposed as 

t~ (res.) t~ (fin-) 

r Vab _ da Vab dcr Vab 



doyab _ da { y^ daf n ' 
dq T dq T dq T 



The first term on the right-hand side contains all the logarithmically-enhanced contributions, 
which have to be resummed to all orders in as, while the second term is free of such contributions 
and can thus be evaluated at fixed order in perturbation theory. Using the Bessel transformation 



between the conjugate variables q? and b (b is the impact parameter), the resummed component 
d&vab can ^ e expressed as 

<i<j (res ' ) M 2 f°° b 

—^(q T ,M, s-a^nDilAilA) = — J db - J {bq T ) W a y 6 (6, M, s; a s (/4),/4/4) , (4) 

where Jq(x) is the Oth-order Bessel function. Considering the Mellin N- moments Wat of W with 
respect to the variable z = M 2 /s at fixed M, the resummation structure of W^ N can be organized 
in exponential form [J 

W£(6, M; a S (/4), a4 &) = H V N (M, a s (/4); M 2 //4 M 2 //4, M 2 /Q 2 ) 

x exp{^(a s (^), L; M 2 //4 M 2 /Q 2 )} , (5) 

were we have defined the logarithmic expansion parameter L = hi(Q 2 b 2 /b"fy, and b Q = 2e~ jE 
{lE — 0.5772... is the Euler number). The scale Q (Q ~ M), which appears on the right-hand side 
of Eq. (jSJ), is the resummation scale [21]. Although W^ (i.e., the product W^ x expj^jy}) does 
not depend on Q when evaluated to all perturbative orders, its explicit dependence on Q appears 
when W^v is computed by truncation of the resummed expression at some level of logarithmic 
accuracy (see Eq. (jBj) below). Variations of Q around M can thus be used to estimate the size of 
yet uncalculated higher-order logarithmic contributions. 

The universale form factor exp{(?7v} contains all the terms that order-by-order in as are loga- 
rithmically divergent as b — > oo (or, equivalently, qr — > 0). The resummed logarithmic expansion 
of the exponent Qn is defined as follows: 

g N (a s , L- M 2 //4, M 2 /Q 2 ) = L g^(a s L) + g {2) (a s L; M 2 //4, M 2 /Q 2 ) 

+ %{JW,M 2 / / 4,M 2 /Q 2 ) + ... (6) 

TV 

where the term Lg^ collects the leading logarithmic (LL) contributions, the function g N includes 
the next-to-leading leading logarithmic (NLL) contributions [17], g N controls the NNLL terms 
[32] |33| [31] and so forth. The explicit form of the functions g^\ g N and g N can be found in 
Ref . [21] . The process dependent function "Hj^ does not depend on the impact parameter b and it 
includes all the perturbative terms that behave as constants as b — > oo. It can thus be expanded 
in powers of «s = ocs(^r)' 



H v N (M,a s ;M 2 /n%M 2 /ti,M 2 /Q 2 ) = a$\M) 1 + ^< (1) (M 2 //4, M 2 /Q 2 ) 

L 7T 

+ (^) 2 < (2) (M 2 /^,M 2 /^,M 2 /Q 2 ) + ...] , (7) 

where <jy is the partonic cross section at the Born level. The first-order coefficients 7~L ^J_ abN 

in Eq. (JTj) are known since a long time [32], while the second-order coefficients Ti q }J. abtN were 
computed only recently [5]. 



'I' For the sake of simplicity we consider here only the case of the diagonal terms in the flavour space of the 
partonic indices a, b. For the general case and a detailed discussion of the resummation formalism, we refer to 
Ref. [21]. 

*The form factor does not depend on the type of produced vector boson. More generally, all the hard-scattering 
processes initiated by quark-antiquark annihilation have the same form factor. 
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Within a straightforward ('naive') implementation of Eq. (EJ), the resummation of the large 
logarithmic contributions would affect not only the small-g^ region, but also the region of large 
values of q?. This can easily be understood by observing that the logarithmic expansion parameter 
L diverges also when b — > 0. To reduce the impact of unjustified higher-order contributions in the 
large-^T region, the logarithmic variable L in Eq. ((5]) is actually replaced by L = In (Q 2 b 2 /bl + 1) 
[211 122] . This replacement has an additional and relevant consequence: after inclusion of the finite 
component (see Eq. (|HJ)), we exactly recover the fixed-order perturbative value of the total cross 
section upon integration of the qx distribution over q? (i.e., the resummed terms give a vanishing 
contribution upon integration over q?). 

We now turn to consider the finite component of the transverse-momentum cross section (see 
Eq. (j3j). Since day does not contain large logarithmic terms in the small-gr region, it can be 
evaluated by truncation of the perturbative series at a given fixed order. In practice, the finite 
component is computed starting from the usual fixed-order perturbative truncation of the partonic 
cross section and subtracting the expansion of the resummed part at the same perturbative order. 
Introducing the subscript f.o. to denote the perturbative truncation of the various terms, we have: 

• (8) 

f.o. 

This matching procedure between resummed and finite contributions guarantees to achieve uni- 
form theoretical accuracy over the region from small to intermediate values of transverse momenta. 
At large values of qx, the resummation (and matching) procedure is eventually superseded by the 
customary fixed-order calculations (their theoretical accuracy in the large-g^ region cannot be 
improved by resummation of the logarithmic terms that dominate in the small-gT region). 

In summary, the inclusion of the functions g^\ g N , H N in the resummed component, to- 
gether with the evaluation of the finite component at LO (i.e. at 0(as)), allows us to perform the 
resummation at NLL+LO accuracy. This is the theoretical accuracy used in our previous study 
[23] of the DY qT distribution. Including also the functions g N and 7i N , together with the 
finite component at NLO (i.e. at 0{a\)) leads to full NNLL+NLO accuracy. The perturbative 
coefficient A^ 3 \ which contributes to the NNLL function g N (see, e.g., Eq. (24) in Ref. [2T]), is 
not yet known. In the following, we assume that the value of A^ is the same as the one [53 ES] 
that appears in resummed calculations of soft-gluon contributions near partonic threshold. Using 
the recently computed H N coefficient [5] , we are thus able to present the complete result for the 
DY ^-distribution up to NNLL+NLO accuracy. We point out that the NNLL+NLO (NLL+LO) 
result includes the full NNLO (NLO) perturbative contribution in the small-gr region. In particu- 
lar, the NNLO (NLO) result for the total cross section is exactly recovered upon integration over 
qx of the differential cross section day /dqx at NNLL+NLO (NLL+LO) accuracy. 

We conclude this section with some comments on the numerical implementation of our calcu- 
lation. Within our formalism, the resummation factor W^(b,M) is directly defined, at fixed M, 
in the space of the conjugate variables b and N. To obtain the hadronic cross section, we have 
to perform inverse integral transformations: the Bessel transformation in Eq. ([3]) and an inverse 
Mellin transformation. These integrals are carried out numerically. The Mellin inversion requires 
the numerical evaluation of some basic iV-moment functions that appear in the expression of the 
the second-order coefficients 'H q } < _ ab N |5J: this evaluation has to be performed for complex values 
of N, and we use the numerical results of Ref. [37J. We recall [21J that the resummed form factor 



exp{G N (as((j, 2 R ),L)} is singular at the values of b where as(^%)L > n/f3 (/3 is the first-order 
coefficient of the QCD /3 function). Performing the Bessel transformation with respect to the 
impact parameter b (see Eq. PJ), we deal with this singularity as we did in Ref. [21] . by using the 
regularization prescription of Refs. [381139]: the singularity is avoided by deforming the integration 
contour in the complex b space. 



3 Numerical results for Z/y* production at the Tevatron 

In this section we consider Z/'-f* production in pp collisions at Tevatron energies. We present 
our resummed results at NNLL+NLO accuracy, we compare them with the NLL+LO results (the 
NLL+LO results in Ref. [21] were obtained by using the MRST2004 NLO parton densities [30]). 
and we comment on the comparison with Tevatron Run II data [301 |3T] . 

The hadronic q? cross section at NNLL+NLO (NLL+LO) accuracy is computed by using the 
MSTW2008 NNLO (NLO) parton densities [H], with a s (/4) evaluated at 3-loop (2- loop) order. 
This choice of the order of the parton densities and as is fully justified both in the small-gr region 
(where the calculation of the partonic cross section includes the complete NNLO (NLO) result 
and is controlled by NNLL (NLL) resummation) and in the intermediate-g^ region (where the 
calculation is constrained by the value of the NNLO (NLO) total cross section). 

As for the electroweak couplings, we use the so called G M scheme, where the input parameters 
are G F , m z , m w . In particular, we use the PDG 2008 [12] values G F = 1.16637 x 10~ 5 GeV -2 , 
mz = 91.1876 GeV, I 'z = 2.4952 GeV, mw = 80.398 GeV. Our calculation implements the decays 
7* — > l + l~ and Z* — > l + l~ at fixed value of the invariant mass of the l + l~ pair. In particular, we 
include the effects of the 7* Z interference and of the finite width of the Z boson. Nonetheless, the 
numerical results presented below are obtained by simply using the narrow-width approximation 
and neglecting the photon contribution. We find that this approximation works to better than 
1% accuracy in the inclusive region of lepton invariant mass that is covered by the DO data_j. 

As discussed in Sect. [21 the resummed calculation depends on the factorization and renormal- 
ization scales and on the resummation scale Q. Our convention to compute factorization and renor- 
malization scale uncertainties is to consider independent variations of fip and /z# by a factor of two 
around the central values Hf = Hr = Tnz (i-e. we consider the range mz/2 < {^f,Hr} < 2mz), 
with the constraint 0.5 < Hf/^r < 2. Similarly, we follow Ref. [21] and we choose Q = mz/2 as 
central value of the resummation scale, considering scale variations in the range m^/4 < Q < mz- 

In Fig. T] (left panel) we present the NLL+LO g-r spectrum at the Tevatron Run II (-^/s = 
1.96 TeVJj. The NLL+LO result (solid line) at the default scales (/i^ = ^r = mz, Q = mz/2) 
is compared with the corresponding LO result (dashed line). The LO finite component of the 
spectrum (see Eq. ([3])), multiplied by a factor of 10 to make it more visible, is also shown for 
comparison (dotted line). We see that the LO result diverges to +00 as qx — > 0. The resummation 
of the small-gr logarithms leads to a well-behaved distribution: it vanishes as qr — > 0, has a 



§The measured qr spectra are inclusive over the following regions of lepton invariant mass: 70-110 GeV [30 
and 65-115 GeV [3T]. 

^Analogous results at the Tevatron Run I (y^i = 1.8 TeV), obtained by using the MRST2004 parton densities 
[40] . were presented in the left panel of Fig. 6 of Ref. 24 . 
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Figure 1: The qx spectrum of Z bosons at the Tevatron Run II: results at NLL+LO (left panel) 
and NNLL+NLO (right panel) accuracy. Each result is compared to the corresponding fixed- order 
result (dashed line) and to the finite component (dotted line) in Eq. |3]). 



kinematical peak at qr ~ 2 GeV, and tends to the corresponding LO result at large values of qT- 
The finite component smoothly vanishes as qx — > and gives a small contribution to the NLL+LO 
result in the low-gr region. 

The results in the right panel of Fig. [I] are analogous to those in the left panel, although 
systematically at one order higher. The qT spectrum at NNLL+NLO accuracy (solid line) is 
compared with the NLO result (dashed line) and with the NLO finite component of the spectrum 
(dotted line). The NLO result diverges to — oo as qx — > and, at small values of q?, it has an 
unphysical peak (the top of the peak is above the vertical scale of the plot) that is produced by 
the numerical compensation of negative leading and positive subleading logarithmic contributions. 
The contribution of the NLO finite component to the NNLL+NLO result is smaller than 1% at 
the peak and becomes more important as qr increases: it is about 8% at qx ~ 20 GeV, about 
20% at q T ~ 30 GeV and about 53% at q T ~ 50 GeV. A similar quantitative behaviour is 
observed by considering the contribution of the NLO finite component to the NLO result; the 
contribution is about 10% at qT ~ 20 GeV, about 22% at qx ~ 30 GeV and about 60% at 
qr ~ 50 GeV. In the region of intermediate values of q? (say, around 20 GeV), the difference 
between the NNLL+NLO and NLO results is larger than the size of the NLO finite component. 
This difference is produced by the logarithmic terms (at NNLO and beyond NNLO) that are 
included in the resummed calculation at NNLL accuracy. At large values of qT the contribution of 
the NLO finite component sizeably increases. This behaviour indicates that the logarithmic terms 
are no longer dominant and that the resummed calculation cannot improve upon the predictivity 
of the fixed-order expansion. 

Comparing the left and right panels of Fig. [I], we see that the qT spectrum is slightly harder 
at NNLL+NLO accuracy than at NLL+LO accuracy. The height of the peak at NNLL+NLO is 
lower than at NLL+LO. This is mainly due to the fact that the NNLO total cross section, which 
fixes the value of the qT integral of our NNLL+NLO result, is only about 3% larger than the NLO 
total cross section, whereas in the region of intermediate values of qT the cross section at NLO 



(and, correspondingly, at NNLL+NLO) is definitely larger than at LO (and, correspondingly, at 
NLL+LO); this leads to a reduction of the cross section at small qx- 
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Figure 2: The qx spectrum of Z bosons at the Tevatron Run II. The bands are obtained by varying 
[ip and fiR (left panel) and Q (right panel) as described in the text. 

In Fig. [2]we show the scale dependence of the NLL+LO (dashed lines) and NNLL+NLO (solid 
lines) results. In the left panel we consider variations of the renormalization and factorization 
scales. The bands are obtained by varying fiR and fip as previously described in this section. 
We note that, in the region of small and intermediate transverse momenta (gT~30 GeV), the 
NNLL+NLO and NLL+LO bands overlap. This feature, which is not present in the case of the 
fixed-order perturbative results at LO and NLO (see Figs. 2 and 3 in Ref. [21]), confirms the 
importance of resummation to achieve a stable perturbative prediction. In the region of small and 
intermediate values of qx, the main difference between the NNLL+NLO and NLL+LO predictions 
is in the size of the scale variation bands. Going from NLL+LO to NNLL+NLO accuracy, we 
observe a reduction of the scale dependence from ±4% to ±3% at the peak, from ±7% to ±3% 
at qx ~ 20 GeV, and from ±7% to ±5% at qr ~ 50 GeV. We point out that the qx region where 
resummed perturbative predictions are definitely significant is a wide region from intermediate 
to relatively-small (say, close to the peak of the distribution) values of qr- In fact, at very small 
values of qx (e.g. gr~5 GeV) the size of non-perturbative effects is expected to be important, 
while in the high-g^ region (e.g. qx ^ 60 GeV) the resummation of the logarithmic terms cannot 
improve the predictivity of the fixed-order perturbative expansion. The inset plot in the left panel 
of Fig. [2] shows the region from intermediate to large values of qx- At large qx, the NLL+LO 
and NNLL+NLO results deviate from each other, and the deviation increases as qx increases. As 
previously stated, this behaviour is not particularly worrying since, in the large-g^ region, the 
resummed results loose their predictivity and can (should) be replaced by customary fixed-order 
results. 

In the right panel of Fig. [2] we consider resummation scale variations. The bands are ob- 
tained by fixing jjlr — ixp — mz and varying Q between m^/4 and mz- Performing variations 
of the resummation scale, we can get further insight on the size of yet uncalculated higher-order 
logarithmic contributions at small and intermediate values of qx- We find that the scale depen- 
dence at NNLL+NLO (NLL+LO) is about ±5% (±12%) in the region of the peak, and about 



7 



±5% (±16%) in the region where q T ~ 20 GeV. We note that in a wide region of qx values, 
5 GeV < qT ^ 50 GeV, the resummation scale dependence is reduced by, roughly, a factor of 2 in 
going from the NLL+LO to the NNLL+NLO result. Comparing the left and right panels of Fig. [2J 
we see that, at NNLL+NLO accuracy, the resummation scale dependence is larger than (though, 
comparable to) the hf and /ir dependence. 

The integral over qT of the resummed NNLL+NLO (NLL+LO) spectrum is in agreement (for 
any values of /ir, /i^ and Q) with the value of the corresponding NNLO (NLO) total cross section 
to better than 1%, thus checking the numerical accuracy of our code. We also note that the 
large-gr region gives a little contribution to the total cross section (see some numerical results in 
Sect. 3.2 of Ref. [24J ) ; therefore, the total cross section constraint mainly acts as a perturbative 
constraint on the resummed spectrum in the region from intermediate to small values of qT- 



0.00 






- 


ii i i i i i i i | i i i i _ 


0,4 


— 


W / NL0 T 


— 


0.2 


— 


\ w, 


T- 




- 


NNLL+NLO * \\ -p 




1' _L - 




~~ J — 4^W 








j 1 ---. 


0.0 
0.2 


x -^K — : vri r 




1 


-L I DO fi + /jT 


_ _ , - -V - 




- 


I DO e + e" \ - 


0.4 


— 


\ — 




- 


: : 1 I : I , 1 1 I 1 I I \ 1" 



10 20 

q T (GeV) 



50 



100 



Figure 3: The normalized qT spectrum of Z bosons at the Tevatron Run II. The NNLL+NLO 
result is compared with the DO data of Refs. l30\ [31\/ . The bands are obtained as described in the 
text. 

The DO Collaboration has measured the normalized qT distribution, ^j 2 -, from data at the 
Tevatron Run II in the e + e~ [30J and yU + /i~ [31] channels. In the left panel of Fig. [3] we report the DO 
data and our corresponding results at NNLL+NLO accuracy. The NNLL+NLO band represents 
our estimate of the perturbative uncertainty, and it is obtained by performing scale variations 
as follows. We independently vary /jlf,/J-r and Q in the ranges niz/2 < {hfiHr} < 2m^ and 
"iz/4 < Q < mz, with the constraints 0.5 < Hf/hr < 2 and 0.5 < Q/fiR < 2. The constraint 
on the ratio {1f/hr is the same as used in the left panel of Fig. [21 it has the purpose of avoiding 
large logarithmic contributions (powers of ln(fi F /fi R )) that arise from the evolution of the parton 
densities. Analogously, the constraint on the ratio Q/hr avoids large logarithmic contributions 
(powers of \n(Q 2 / ft, R )) in the perturbative expansion of the resummed form factoiUJI exp{Q N } (see 
Eq. (jSJ)). We recall (see e.g. Eq. (19) of Ref. [21]) that the exponent Qn of the form factor is 
obtained by q 2 integration of perturbative functions of cio,(q 2 ) over the range b^/b 2 < q 2 < Q 2 . 
To perform the integration with systematic logarithmic accuracy, the running coupling cts(q 2 ) is 
then expressed in terms of ag(/z#) (and ln(g 2 //i^)). As a consequence, the renormalization scale 

"We do not apply additional constraints on the ratio Q/ /ip, since the form factor does not depend on fip. 



fj,R should not be too different from the resummation scale Q, which controls the upper bound of 
the q 2 integration. 

The DO data and the NNLL+NLO band are presented in the left panel of Fig. [3j The inset 
plot shows the region from q T = 30 GeV up to qx = 100 GeV. A quick inspection of the figure 
shows that the data are described quite well by the NNLL+NLO perturbative predictions. 

Differences and similarities between theoretical calculations and the data are more clearly 
visible by considering their fractional difference with respect to a 'reference' theoretical result. We 
choose the NNLL+NLO result at central values of the scales (i.e. \ip = Hr = rriz, Q = rnz/2) 
as 'reference' theory, and we show the ratio (X— theory) /theory in the right panel of Fig. [31 The 
label X refers to either the experimental data or the NNLL+NLO (solid lines) and NLO (dashed 
lines and dotted line) results, including their scale dependence. 

Considering the right panel of Fig. [3, we first comment on the scale uncertainty band of the 
NNLL+NLO result (solid lines). Such uncertainty is about ±6% at the peak, it decreases to about 
±4-5% in the region up to qx = 10 GeV, and then it increases, reaching the size of about ±12% at 
qT = 50 GeV. In the region beyond qx ~ 60 GeV the resummed result looses predictivity, and its 
perturbative uncertainty becomes large. The right panel of Fig. [3] also shows the scale variation 
band of the NLO result. The NLO band (dashed lines) is obtained by varying /j, F and hr (the NLO 
calculation does not depend on the resummation scale Q) as in the NNLL+NLO calculation. We 
comment on the comparison between the NLO and NNLL+NLO bands. At large valueo of qr, 
the NLO and NNLL+NLO bands overlap (the NLO and NNLL+NLO are certainly consistent), 
and the NLO result has a smaller uncertainty. At intermediate values of transverse momenta, the 
NLO result is lower than the NNLL+NLO result, and the corresponding scale variation bands 
do not overlap. We recall (see the discussion in Sect. 3.1 of Ref. [21]) that in this region the 
NLO band underestimates the true perturbative uncertainty of the NLO result; indeed, the NLO 
band and the corresponding LO band do not overlap at intermediate values of q F (see Figs. 3 and 
5 in Ref. [23]). To get some quantitative insight into the 'true' perturbative uncertainty of the 
NLO calculation in this q F region, we can consider wider scale variations and, in particular, we 
can lower the values of fip an d fiR. In the right panel of Fig. [31 we show the NLO band that 
we obtain by independently varying fip and {ir in the range m^/4 < {iif.^r} < 2m^, with the 
constraint 0.5 < Hf/hr < 2: this band is delimited by the dotted line and the lower dashed line 
(the region between the dotted line and the central values of the dashed band roughly corresponds 
to scale variations in the range m z /4 < {fj, F ,p, R } < m z)- We note that lowering the scales at 
NLO improves the consistency between the NLO and NNLL+NLO results. We also note that we 
have considered similar enlarged scale variations at NNLL+NLO accuracy, and we have checked 
that they do not significantly modify the NNLL+NLO band in the intermediate-gT region. This 
confirms the better stability of the NNLL+NLO calculation with respect to scale variations. In 
the small-gr region, the NLO result is theoretically unreliable. The NLO band quickly deviates 
from the NNLL+NLO band as q F decreases. 

The right panel of Fig. [3] shows that the NNLL+NLO result is consistent with the DO data, 
even at very low values of q F - We note that the resummed result is obtained in a perturbative 
framework. At low values of qr, non-perturbative effects are important and are expected (see, 
e.g., the final part of Sect. 5 in Ref. [23]) to shift the resummed result such as to improve the 



**The available DO data in the region 100 GeV< qr < 250 GeV are consistent with the NLO result (see 
Refs. [301 E] and Fig. 5 in Ref. [M]). 
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agreement with the data. In the region where qx ^ 50 GeV, the experimental errors and the 
corresponding NNLL+NLO errors overlap, with the sole exception of a couple of data points at 
very low q?. In the same region, the perturbative uncertainty of the NNLL+NLO result turns out 
to be comparable with the size of the experimental errors. As pointed out by the DO Collaboration 
[30| [31] . the NLO result tends to undershoot the data in the region of intermediate values of qr'- 
NNLL resummation improves the agreement with the data in this q? region. 



4 Summary 

In this paper we have considered the qx spectrum of DY lepton pairs produced in hadron collisions, 
and we have presented a perturbative QCD study based on transverse-momentum resummation 
at the NNLL order. 



We have followed the formalism developed in Refs. [201 121] , which is valid for the production of a 
generic high-mass system of non strongly-interacting particles in hadron collisions. The formalism 
combines small-g^ resummation at a given logarithmic accuracy with the fixed-order calculations. 
It implements a unitarity constraint that guarantees that the integral over q? of the differential 
cross section coincides with the total cross section at the corresponding fixed-order accuracy. This 
leads to QCD predictions with a controllable and uniform perturbative accuracy over the region 
from small up to large values of qx- At large values of qr, the resummation formalism is superseded 
by customary fixed-order calculations. 

We have considered the explicit case of DY lepton pairs from the decay of a Z boson produced at 
the Tevatron Run II. Using the recently computed NNLL coefficient 7i N [5], we have extended 
the NLL+LO resummed calculation presented in Ref. [21] to the NNLL+NLO accuracy. The 
NNLL corrections are not large and make the gr spectrum slightly harder. We have performed 
a study of the scale dependence of the calculation to estimate the corresponding perturbative 
uncertainty. In a wide region of transverse momenta (5 GeV < qr ^ 50 GeV) the size of the scale 
uncertainties is considerably reduced in going from NLL+LO to NNLL+NLO accuracy. 



We have compared the resummed calculation with the results of measurements [3UI EI] of the 
normalized q? spectrum at the Tevatron Run II. The perturbative uncertainty of the NNLL+NLO 
results turns out to be comparable with the experimental errors. The NNLL+NLO results (without 
the inclusion of any non-perturbative effects) are consistent with the experimental data in a wide 
region of transverse momenta. Comparing the NNLL+NLO and NLO results, we have also shown 
that NNLL resummation improves the agreement with the data at intermediate values of gr- As is 
well known (and theoretically expected), the NLO result fails to describe the data at small values 
of qx- 

More detailed comparisons with available data on vector boson production and further studies 
of theoretical uncertainties, including the impact of non-perturbative effects, are left to future 
investigations. 

Note added. After the completion of this paper, the value of the coefficient A^ (see the related 
comment in the final part of Sect. [2]) for qx resummation was derived in Ref. [13]. We have checked 
the quantitative effect of this value of A^ on our results for the qx distribution of Z bosons. We 
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find that the effect is generally very small. The largest effect is produced in the region of very low 
values of qx ; for instance, in the case of Z production at the Tevatron, the quantitative effect is 
at the level of about 2% (4%) at q T ~ 2 GeV (1 GeV). 
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